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ABSTRACT 


Isentropic  potential  vorticity  (IPV)  fields  calculated  from  the  Navy  Operational 
Regional  Atmospheric  Prediction  System  are  analyzed  to  determine  their  usefulness  as 
an  aid  to  tropical  cyclone  recurvature  forecasts.  The  IPV  fields  associated  with  Typhoon 
Nelson  are  calculated  on  a  80  km  grid  for  the  period  0000  UTC  4  October  to  0000  UTC 
7  October  1988.  In  this  preliminary  study,  IPV  advection  fields  and  the  horizontal  and 
vertical  structure  of  the  bogus  vortex  are  examined  to  determine  their  contributions  to 
the  model  forecast.  The  extremely  broad  bogus  vortex  is  found  to  distort  the  IPV  fields 
and  cause  regions  of  negative  IPV,  which  is  believed  to  result  in  inertial  instability. 
Approximately  36  h  are  required  for  the  model  to  adjust  to  the  initial  imbalance  in  the 
mass  and  wind  fields.  It  appears  that  the  unrealistic  bogus  vortex  representation  may 
have  degraded  the  early  portions  of  the  model  forecast.  Analysis  of  IPV  advection  fields 
each  12  h  during  the  NORAPS  model  integration  showed  that  storm  movement  is  pri¬ 
marily  due  to  the  influence  of  self-advection  and  the  large-scale  steering  flow.  During 
later  periods  other  advective  features  associated  with  adjacent  synoptic  systems  begin  to 
influence  the  storm  movement.  Forecasters  with  access  to  such  IPV  fields  may  be  able 
to  evaluate  the  likely  validity  of  the  model  forecast  of  recurvature.  The  usefulness  of  the 
IPV  representation  from  the  NORAPS  prediction  in  this  single  case  study  suggest  that 
additional  cases  of  tropical  cyclone  recurvature  be  examined  in  terms  of  IPV  concepts. 
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I.  INTRODUCTION 


Tropical  cyclones,  which  account  for  the  strongest  sustained  winds  observed,  pose 
an  extreme  hazard  to  military  and  civilian  operations.  Inaccurate  forecasts  of  the 
movement  of  these  storms  have  led  to  loss  of  life  and  resources.  According  to  Elsberrv 
(1987),  tropical  cyclone  forecasts  have  improved  relative  to  the  1960's  in  terms  of  re¬ 
ducing  error  magnitude  as  well  as  in  terms  of  consistency.  Long-term  forecasts  (48  to 
72  h)  have  the  greatest  errors,  especially  for  storms  that  have  erratic  tracks  and  for  those 
that  recurve.  Elsberrv  (1987)  states  that  track  errors  associated  with  an  incorrect  as¬ 
sessment  of  the  recurvature  situation  may  exceed  1850  km  (1000  n  mi)  at  72  h.  Because 
the  impact  of  these  forecasts  is  usually  very  high,  it  is  essential  that  techniques  be  de¬ 
veloped  to  improve  recurvature  forecasts. 

An  accurate  long-term  forecast  is  critical  for  military  operations.  The  time  required 
to  secure  or  deploy  resources  to  prevent  damage  due  to  high  winds  is  beyond  the  time 
period  when  the  most  success  has  been  achieved  in  tropical  cyclone  forecasting.  An  in¬ 
accurate  long-term  forecast  can  have  very  costly  consequences.  Loss  of  life,  aircraft  and 
mission  readiness  may  result  from  these  deficiencies. 

The  goal  of  this  study  is  to  demonstrate  the  utility  of  potential  vorticity  fields  during 
the  recurvature  of  tropical  cyclones  as  they  interact  with  midlatitude  cyclones. 
Isentropic  potential  vorticity  (IPV)  fields  associated  with  both  the  tropical  and 
midlatitude  cyclone  are  analyzed  to  determine  if  the  orientation  and  advection  of  these 
fields  with  respect  to  the  two  systems  may  be  a  precursor  for  the  recurvature  of  the 
tropical  cyclone. 

IPV  fields  have  provided  a  useful  description  of  extratropical  cyclogenesis  (Hoskins 
et  al.  1985).  Other  researchers  (Uccellini  el  al.  1985)  have  made  use  of  IPV  fields  to  de¬ 
scribe  explosive  cyclogenesis.  Normally,  IPV  would  not  be  considered  relevant  in  the 
quasi-barotropic  tropics.  However,  the  important  interactions  with  extratropical  sys¬ 
tems  during  recurvature  may  be  described  well  in  terms  of  IPV  lobes  approaching  the 
tropical  cyclone  such  that  the  advection  of  IPV  accounts  for  the  deflection  of  the  trop¬ 
ical  cyclone  into  the  midlatitudes. 

Just  as  the  sparsity  of  observations  has  hindered  studies  of  maritime  extratropical 
cyclones,  the  lack  of  observations  also  impedes  the  use  of  IPV  in  the  tropics.  Avail¬ 
ability  of  data  only  at  the  mandatory  pressure  levels,  and  thus  only  poor  vertical  resol- 


ution,  has  made  it  extremely  difficult  to  describe  accurately  the  vertical  structure  oflPV 
in  tropical  cyclones.  Researchers  have  overcome  these  hurdles  by  utilizing  various 
models.  Schubert  and  Alworth  (1987)  and  Thorpe  (1985)  make  use  of  analytical  models 
to  describe  IPV  structures  associated  with  the  tropical  cyclone  vortex.  The  use  of  nu¬ 
merical  models  in  determining  the  vertical  structure  of  IPV  in  extratropical  cyclones  has 
been  handicapped  due  to  the  lack  of  vertical  resolution  (Elsberry  and  Kirchoffer  1988). 

This  study  describes  IPV  features  using  a  a  research  version  of  the  Navy  Operational 
Regional  Atmospheric  Prediction  System  (NORAPS).  NORAPS  has  an  advantage  over 
other  models  in  that  its  horizontal  (80  km)  and  vertical  (21  levels)  resolution  should  be 
sufficient  to  describe  accurately  the  small-scale  features.  Thus,  an  analysis  of  the 
NORAPS  data  should  show  interaction  between  the  midlatitude  and  tropical  cyclone. 
Typhoon  Nelson,  the  only  super-typhoon  of  the  19S8  season,  was  selected  because  it 
was  accurately  forecast  by  NORAPS.  The  results  of  this  study  may  facilitate  the  use 
of  operational  IPV  analysis  and  hopefully  allow  the  forecaster  to  evaluate  the  4S  and  72 
h  forecasts  of  recurving  tropical  cyclones. 

Chapter  II  summarizes  past  research  with  IPV  fields.  Chapter  III  describes  the 
synoptic  situation  during  the  recurvature  of  Nelson  and  also  summarizes  the  NORAPS 
forecast  of  this  storm.  A  description  of  Typhoon  Nelson  in  terms  of  IPV  fields  is  pre¬ 
sented  in  Chapter  IV.  This  chapter  also  includes  a  section  on  the  advection  of  potential 
vorticity  and  how  this  may  affect  the  propagation  of  the  storm.  Appendix  A  gives  an 
indepth  description  of  the  NORAPS  model,  and  Appendix  B  outlines  the  procedures  for 
the  calculation  of  potential  vorticity  from  the  NORAPS  fields.  Conclusions  and  sug¬ 
gestions  for  further  study  are  outlined  in  Chapter  V. 
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II.  BACKGROUND 


The  significance  of  potential  vorticity  (PV)  and  its  effect  on  meteorological  phe¬ 
nomena  cannot  be  over-emphasized.  Within  the  last  10  years,  researchers  have  at¬ 
tempted  to  explain  such  phenomena  as  cyclogenesis  and  explosive  cyclogenesis  in  terms 
of  the  distribution  and  structure  of  PV  in  the  atmosphere.  Hoskins  et  al.  (1985)  provide 
a  historical  perspective  for  a  number  of  early  studies  concerning  the  details  and  upper-air 
structure  of  PV. 

1.  Midlatitude  applications 

The  conservative  properties  of  PV  make  it  an  ideal  tracer  of  air  parcels  in  the 
atmosphere.  Reed  (1955)  and  Reed  and  Sanders  (1953)  made  use  of  this  concept  to 
distinguish  air  that  is  of  stratospheric  origin  from  tropospheric  air  due  to  the  large  values 
of  PV  normally  found  in  the  stratosphere.  The  greater  thermal  stability  of  the 
stratosphere  provides  for  the  high  values  of  PV.  Shapiro  (19S0)  confirms  the  application 
of  this  concept  in  his  calculations  of  PV  on  isentropic  surfaces  from  aircraft  observations 
of  jet  streams. 

Reed  (1955)  applied  the  PV  tracer  by  demonstrating  that  stratospheric  air  could 
be  advected  down  sloping  isentropic  surfaces  to  altitudes  well  within  the  troposphere. 
This  injection  of  stratospheric  air  into  the  tropopause  is  defined  as  a  tropopause  fold. 
This  phenomenon  is  a  mechanism  for  detaching  stratospheric  air  with  high  values  of  PV 
and  transporting  this  air  toward  the  mid-troposphere  along  the  axis  of  the  jet.  Reed 
suggested  that  there  is  a  link  between  the  tropopause  fold  and  upper-level  frontogenesis. 
Daniclsen  (196S)  confirmed  Reed's  hypothesis  and  stated  that  the  large  values  of  PV 
normally  observed  in  upper-level  fronts  are  of  stratospheric  origin.  Tropopause  folds 
usually  occur  to  the  west  of  developing  troughs  along  the  axis  of  confluence  separating 
the  thermodynamic  direct  and  indirect  circulation  cells  beneath  the  core  of  the  jet 
(Danielsen  1968).  Upper-level  fronts  that  are  located  in  the  negative  vorticity  region 
between  the  trough  and  the  upstream  ridge  are  located  in  a  region  of  downward  motion 
where  the  strongest  subsidence  occurs  near  the  warm  edge  of  the  front.  The  sinking  of 
warm  air  intensifies  the  temperature  gradient  due  to  adiabatic  compression.  In  his  study 
of  a  midlatitude  cyclone,  Shapiro  (1980)  concluded  that  intense  baroclinic  processes 
within  and  on  the  scale  of  upper-level  fronts  are  a  key  factor  in  extratropical 
cyclogenesis. 
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Uccellini  et  al.  (1985)  propose  that  the  high  values  of  isentropic  potential 
vorticity  (IPV)  found  in  the  tropopause  fold  and  associated  with  the  upper-level  front 
may  be  a  possible  triggering  mechanism  for  the  explosive  development  of  extratropical 
cyclones.  They  suggest  that  cyclogenesis  could  occur  as  stratospheric  air  with  high  val¬ 
ues  of  PV  is  injected  into  the  troposphere,  and  the  vortex  tubes  are  then  stretched  as  the 
static  stability  decreases  significantly.  If  the  high  values  of  stratospheric  PV  are  to  be 
conserved  then  the  absolute  vorticity  must  increase,  which  will  enhance  the  cyclonic 
circulation  in  the  troposphere.  If  this  upper-level  vorticity  maximum  is  superposed  on 
a  region  of  low-level  baroclinity,  cyclogenesis  may  occur. 

Uccellini  et  al.  (1985)  demonstrate  that  the  injection  of  high  values  of  PV  into 
the  lower  levels  of  the  troposphere  may  have  been  a  triggering  mechanism  for  the  ex¬ 
plosive  development  of  Presidents'  Day  cyclone  of  18-19  February  1979.  They  show  that 
the  high  values  of  PV  were  advected  over  the  storm  center  when  explosive  cyclogenesis 
occurred.  They  concluded  that  the  explosive  development  of  the  cyclone  was  enhanced 
by  the  descent  of  stratospheric  air  and  the  extension  of  high  potential  vorticity  to  lower 
levels. 

Although  the  usefulness  of  IPV  fields  for  understanding  the  structure  and  be¬ 
havior  of  mid-latitude  weather  systems  has  been  demonstrated,  operational  use  of  IPV 
fields  has  been  hampered  by  various  factors  (Elsberrv  and  Kirchoffer  198S).  Before  the 
advent  of  computers,  the  extra  computational  time  needed  to  create  isentropic  analysis 
relative  to  isobaric  analysis  was  a  hindrance.  The  second  factor  is  the  need  for  signif¬ 
icant  level  data  to  accurately  define  the  positions  and  magnitudes  associated  with  fea¬ 
tures  such  as  tropopause  folds,  inversions  and  other  stable  layers.  This  restriction  is 
critical  for  maritime  analysis  where  few  rawinsondes  exist.  The  lack  of  observations 
hinders  the  use  of  the  IPV  tools  in  conjunction  with  tropical  cyclone  structure.  To 
overcome  this  inhibiting  factor,  researchers  have  used  numerical  and  analytical  models. 
Operational  numerical  models  do  not  have  adequate  vertical  resolution  to  describe  key 
IPV  features.  Numerical  models  may  not  realistically  represent  all  of  the  physical 
processes  involved. 

2.  Tropical  cyclone  application 

Various  analvtical  models  have  been  used  to  simulate  the  vortex  structure  of 
tropical  cyclones  in  terms  of  PV.  Thorpe  (1985)  uses  PV  concepts  in  extending  quasi- 
geostrophy  and  semi-geostrophy  to  balanced  curved  flows.  His  equation  set  consists  of 
a  prognostic  equation  for  the  PV  and  two  diagnostic  equations  for  the  potential  function 
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where  </>  is  the  geopotential  and  v  is  the  gradient  wind.  These  equations  give  the  bal¬ 
anced  vortex  structure  and  an  equation  for  the  agradient  circulation  produced  by 
diabatic  forcing.  Thorpe  (1985)  indicates  that  the  vortex  structure  is  completely  deter¬ 
mined  at  any  given  time  by  the  solution  of  the  potential  function  equation,  if  the 
boundary  variations  of  potential  temperature  and  the  interior  variation  of  PV  are 
known.  He  concludes  that  the  variations  in  PV  can  occur  either  in  association  with  a 
sloping  tropopause  or  in  the  troposphere  due  to  diabatic  heating  that  produces  both 
sources  and  sinks  of  PV.  A  tropopause  that  descends  over  the  vortex  produces  a  nega¬ 
tive  radial  gradient  of  PV  that  is  consistent  with  cyclonic  flow.  This  phenomena  would 
counteract  the  tendency  for  upper  anticyclonic  circulation  due  to  the  warm  core  struc¬ 
ture  of  the  tropical  cyclone.  He  determined  that  the  inclusion  of  a  tropopause  in  the 
model  also  increased  the  low-level  cyclonic  flow.  Conversely,  a  tropopause  that  is  ele¬ 
vated  in  the  core  is  consistent  with  anticyclonic  flow.  Thorpe's  first  attempts  to  describe 
PV  structure  in  a  tropical  cyclone  using  an  analytical  model  proved  to  be  the  starting 
point  for  other  researchers. 

Schubert  and  Alworth  (1987)  also  make  use  of  an  analytical  model  in  their  de¬ 
scription  of  the  evolution  of  PV  in  tropical  cyclones.  They  solve  Eliassen's  (1952)  bal¬ 
anced  vortex  equations  analytically  so  that  the  value  of  PV  is  known  at  all  times. 
Making  use  of  the  "invertibility  principle"  as  described  by  Hoskins  et  al.  (1985),  they  di¬ 
agnostically  deduce  the  complete  flow  structure  of  the  tropical  cyclone  from  the  spatial 
distribution  of  PV.  Their  numerical  results  indicate  that  latent  heat  release  produces  a 
region  of  large  PV  in  lower  levels  and  a  region  of  small  PV  at  upper  levels.  Due  to  ver¬ 
tical  advection,  the  upper  tropospheric  minimum  is  pinched  off  as  large  PV  spreads 
throughout  the  troposphere.  This  configuration  causes  tangential  winds  to  develop  at 
an  increasing  rate,  so  that  typhoon  strength  winds  are  attained  within  96  h.  If  the  PV 
is  allowed  to  continue  to  evolve,  intense  low-level  winds  are  produced. 

Schubert  and  Alworth  (19S7)  and  Thorpe  (1985)  introduce  a  somewhat  different 
view  of  tropical  cyclones  with  primary'  attention  focused  on  the  PV  field.  Since  PV  fields 
are  very  significant  in  midlatitude  baroclinic  waves  and  frontogenesis,  researchers  have 
moved  closer  to  the  goal  of  unified  dynamical  framework  for  the  study  of  midlatitude 
and  tropical  cyclones. 
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The  lack  of  high-resolution  observations  is  a  great  problem  in  determining  the 
structure  of  IPV  in  tropical  cyclones.  Earlier  studies  used  analytical  models  to  overcome 
these  resolution  problems.  This  study  attempts  to  represent  the  IPV  structure  in  a 
tropical  cyclone  and  a  midlatitude  cyclone  during  recurvature  by  using  numerical  model 
output.  The  Navy  Operational  Regional  Atmospheric  Prediction  System  (NORAPS) 
output  fields  are  selected  to  describe  the  key  upper-level  features  because  of  its  increased 
vertical  resolution. 

Comparisons  of  key  IPV  structures  is  the  goal  of  this  study.  As  mentioned  by 
Hoskins  etal.  (1985),  IPV  maxima  associated  with  stratospheric  air  in  the  tropopause 
fold  is  a  key  feature  in  the  upper-level  cyclone  whereas  a  maximum  (minimum)  of  IPV 
is  expected  in  the  low  levels  (upper  levels)  of  the  tropical  cyclone.  Therefore,  it  is  ex¬ 
pected  that  an  interaction  due  to  the  IPV  gradients  should  occur  as  the  tropical  cyclone 
moves  poleward.  This  study  attempts  to  show  that  the  interaction  of  the  two  lobes  of 
IPV  and  how  they  are  oriented  in  relation  to  the  two  systems  may  be  a  precursor  in 
forecasting  the  recurvature  of  a  tropical  cyclone. 


6 


III.  SYNOPTIC  DISCUSSION  AND  FORECAST  TRACK  OF  TYPHOON 

NELSON 


A.  SYNOPTIC  DISCUSSION 

Typhoon  Nelson  was  the  first  major  tropical  cyclone  in  the  western  North  Pacific 
during  October  1988.  A  lull  in  tropical  cyclone  activity  had  occurred  in  the  later  part 
of  September  as  polar  air  pushed  southward  across  the  Asian  mainland  and  Japanese 
Islands.  During  this  period,  the  monsoon  trough  had  returned  to  its  normal  position. 
Nelson  developed  in  the  monsoon  trough  over  the  Philippine  Sea  and  initially  moved 
toward  the  west  and  then  northwest  as  it  tracked  along  the  periphery  of  the  subtropical 
ridge.  Nelson  intensified  as  it  moved  to  the  northwest  and  reached  peak  winds  of  140 
kt  (72  m  s)  at  1200  UTC  4  October.  By  0000  L'TC  5  October,  Nelson  had  begun  to  re¬ 
curve  (Fig.  1).  As  will  be  described  below,  the  synoptic  situation  during  this  time  period 
was  ideal  for  recurvature  to  the  northeast.  An  upper-level  trough  located  between  the 
stationary  Tibetan  High  and  the  subtropical  high  over  the  western  Pacific  deflected 
Nelson  to  the  northeast.  The  relatively  stationary  large-scale  features  controlling  the 
motion  probably  contributed  to  an  accurate  forecast  of  Nelson's  recurvature  (to  be  dis¬ 
cussed  below). 

Key  aspects  of  the  synoptic  evolution  that  relate  to  the  movement  and  the  inter¬ 
action  with  a  midlatitude  cyclone  are  highlighted  in  the  next  section.  Two  levels  are  se¬ 
lected  for  analysis.  The  upper  level  (300  mb)  provides  an  accurate  description  of  the 
interaction  of  Nelson  with  the  midlatitude  cyclone,  whereas  the  midlevel  (500  mb)  is  se¬ 
lected  to  represent  the  steering  flow.  Although  considerable  debate  still  exists  as  to 
which  level  or  layer  most  accurately  represents  the  tropical  cyclone  movement,  a  high 
correlation  exists  between  tropical  cyclone  displacement  and  the  flow  at  the  mid- 
tropospheric  levels.  Elsberry  (1987)  states  that  the  best  representation  may  be  a  deep- 
layer  mean  flow  that  would  include  winds  of  850,  700,  500,  400  and  250  mb.  In  this  case, 
an  adequate  description  can  be  achieved  with  the  500  mb  winds. 

NOGAPS  analysis  in  the  immediate  region  of  the  typhoon  is  suspect  due  to  the  lack 
of  observations.  Key  typhoon  features,  such  as  upper-level  cyclonic  flow  near  the  center 
and  a  transition  to  anticyclonic  flow  at  larger  radius,  are  not  described  and  must  be 
bogused  into  the  analysis.  The  analysis  is  also  limited  in  horizontal  resolution.  Small- 
scale  features  occurring  during  the  interaction  of  the  typhoon  with  the  midlatitude 
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cyclone  are  not  represented  well  due  to  the  sparsity  of  observations.  To  overcome  this 
limitation,  \ORAPS  data  with  a  resolution  of  80  km  will  be  used  later  to  describe  these 
features. 

1.  0000  UTC  4  October  1988 

The  post-storm  analysis  (best  track)  has  Typhoon  Nelson  located  at  18.5°N, 
125. 8°E  with  winds  of  120  kt.  The  storm  track  is  to  the  northwest  at  a  speed  of  11-12 
kt.  Significant  300  mb  features  during  this  period  (Fig.  2a)  include  a  trough  that  extends 
from  31°N,  120°E  to  24°N,  119°E.  The  position  of  the  trough  with  respect  to  Nelson 
is  advantageous  in  the  rapid  intensification  or  sustained  development  of  the  tropical 
cyclone.  According  to  Elsberry  (1987),  the  development  of  a  strong  outflow  channel  to 
the  westerlies  as  an  upper-level  trough  approaches  or  develops  to  the  west  of  the  cyclone 
is  advantageous  for  system  development.  The  outflow  jet  provides  a  good  exhaust  sys¬ 
tem  for  the  storm  and  thus  enhances  intensification.  A  col  area  is  located  just  west  of 
Taiwan  between  two  anticyclones:  one  positioned  east  of  the  storm  at  21°N,  136.5°E, 
and  the  other  located  over  southern  China.  The  anticyclone  east  of  the  storm  provides 
good  outflow,  so  continued  intensification  is  expected. 

The  midlevel  analysis  (Fig.  2b)  is  dominated  by  a  high  pressure  center  located 
at  23.5°N,  113°E  and  the  subtropical  ridge  (STR)  located  over  the  western  Pacific.  The 
500  mb  flow  indicates  movement  along  the  periphery ’  of  the  STR.  However,  the  move¬ 
ment  and  intensity  of  the  high  pressure  cell  over  southern  China  is  also  a  factor  deter¬ 
mining  the  direction  in  which  Nelson  will  propagate.  Typhoon  Nelson  appears  to  be 
moving  toward  the  break  in  the  ridge  northwest  of  the  storm  postion,  which  is  repres¬ 
ented  by  a  col  area  northeast  of  Taiwan. 

2.  1200  UTC  4  October  1988 

Typhoon  Nelson  is  now  located  at  20°N,  124. 7°E  and  has  peak  winds  of  140 
kt.  Nelson  is  tracking  through  an  area  where  a  large  number  of  tropical  cyclones  reach 
supertyphoon  intensity,  according  to  the  1970  Joint  Typhoon  Warning  Center.  Al¬ 
though  the  storm  continues  to  track  to  the  northwest,  its  speed  of  propagation  has  de¬ 
creased  to  7-8  kt. 

A  key  300  mb  feature  is  the  location  of  a  weakness  in  the  ridge  line  with  respect 
to  Nelson  (Fig.  3a).  This  feature  extends  from  30°N,  121°E  to  25°N,  121°E.  The  strong 
winds  in  excess  of  50  kt  to  the  northwest  of  the  storm  center  provide  a  good  outflow  jet, 
which  aids  in  the  intensification  of  the  tropical  cyclone.  The  upper-level  col  located  at 
24°N,  1 2 1  °E  has  propagated  very  slowly  to  the  east.  However,  precise  displacement  of 
this  feature  is  in  doubt  due  to  the  sparsity  of  observations  in  this  area. 
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Figure  2.  NOGAPS  analysis  for  00  UTC  4  Oct  1988  at(A)  300  mb  and  (B)  500 
mb.  Streamlines  are  analyzed  in  the  region  of  Typhoon  Nelson  (dot)  and  the  sub¬ 
tropical  highs. 


An  analysis  of  the  mid-level  flow  provides  the  reasoning  for  the  slow  propa¬ 
gation  of  Typhoon  Nelson  (Fig.  3b).  Nelson  continues  to  track  around  the  periphery 
of  the  STR  towards  the  area  of  least  resistance.  The  storm  is  located  in  an  area  of  weak 
mid-level  winds  just  south  of  the  col  between  the  high  cell  located  at  22°N,  Ill  °E  and 
the  STR.  The  northeast-southwest  orientation  of  the  subtropical  ridge  axis  contributes 
to  more  meridional  steering  flow. 

3.  0000  UTC  5  October  1988 

During  this  period,  Typhoon  Nelson  has  begun  to  recurve  to  the  northeast.  The 
storm  is  located  at  21.3°N,  124.4°E  with  winds  of  140  kt.  The  300  mb  circulation  shows 
slight  troughing  to  the  northwest  of  the  storm  (Fig.  4a).  Strong  winds  on  the  down¬ 
stream  side  of  the  trough  and  on  the  periphery  of  the  western  Pacific  ridge  continue  to 
be  in  the  proper  position  to  aid  in  the  intensification  of  Nelson.  The  col  located  over 
Taiwan  and  the  anticyclone  located  over  southern  China  have  remained  relatively  sta¬ 
tionary  in  the  last  12  h. 

The  movement  of  Nelson  has  slowed  to  5  kt  as  recurvature  is  occurring.  The 
500  mb  analysis  shows  that  the  high  over  southern  China  is  building  (Fig.  4b).  The 
winds  between  the  two  highs  in  the  region  of  the  storm  are  opposing,  which  may  be  a 
factor  in  the  slow  propagation  of  Nelson.  The  col  northeast  of  Taiwan  has  remained 
stationary  as  Nelson  continues  to  track  toward  the  area  of  least  resistance. 

4.  1200  UTC  5  October  1988 

Typhoon  Nelson  is  located  at  22.5°N,  125.1°E  with  winds  of  140  kt.  Nelson  is 
propagating  to  the  northeast  at  7-8  kt.  The  upper-level  analysis  shows  a  midlatitude 
trough  located  at  30°N,  115°E  is  approaching  from  the  northwest  (Fig.  5a).  The 
strongest  winds  associated  with  the  downstream  side  of  the  trough  are  located  to  the 
northwest  of  the  storm  center.  This  flow,  in  addition  to  the  winds  associated  with  the 
western  Pacific  ridge,  still  provide  an  excellent  upper-level  outflow.  The  col  located  over 
Taiwan  has  remained  stationary  in  the  last  12  h.  Although  the  trough  and  Nelson  are 
quite  close,  Nelson  continues  to  be  a  warm  core  system  with  little  or  no  penetration  of 
cold  air  from  the  north. 

Analysis  of  the  mid-level  features  (Fig.  5b)  shows  a  strong  tendency  for  con¬ 
tinued  movement  toward  the  northeast.  The  circulation  pattern  has  not  changed  sig¬ 
nificantly  in  the  last  12  h.  The  steering  flow  near  the  storm  center  is  to  the  northeast 
along  the  periphery  of  the  STR.  The  location  of  the  col  and  the  high  over  southern 
China  are  unchanged.  The  only  significant  change  is  a  slight  intensification  of  the  STR, 
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Figure  3.  As  in  Fig.  2,  except  for  12  UTC  4  Oct  1988. 

although  the  analyzed  change  in  intensity  may  be  suspect  due  to  the  sparsity  of  obser¬ 
vations  in  this  region. 
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5.  0000  UTC  6  October  1988 

Typhoon  Nelson  has  continued  to  track  to  the  northeast  at  10  kt,  and  it  is  now 
located  at  24.0°N,  126.7°E  and  has  weakened  with  peak  winds  of  125  kt.  The  upper-level 
analysis  (Fig.  6a)  shows  cyclonic  circulation  has  developed  around  Nelson.  The  col  area 
now  located  at  27°N,  126°E  has  slowly  drifted  to  the  east  and  is  located  directly  north 
of  the  storm.  The  trough  has  also  propagated  to  the  east  as  it  extends  from  32°N, 
126°E  to  27°N,  126.2°E.  The  trough  remains  west  of  the  cyclone  and  is  becoming  more 
north-south  aligned.  Although  the  typhoon  is  getting  closer  to  the  midlatitude  trough, 
the  systems  remain  separate.  The  downstream  side  of  the  trough  still  provides  a  good 
outflow.  However,  the  upper-level  winds  in  this  region  have  weakened  in  the  last  12  h 
due  to  the  orientation  of  the  trough  with  respect  to  the  position  of  the  tropical  cyclone. 
Nelson  begins  to  lose  the  upper-level  support  that  is  so  critical  for  storm  development 
and  maintenance. 

The  mid-level  analysis  (Fig.  6b)  shows  the  midlatitude  trough  has  deepened  as 
it  extends  from  34°N,  131°E  to  28°N,  128°E.  The  col  area  has  propagated  to  the  east 
and  is  located  directly  north  of  the  cyclone.  Strong  (50  kt)  winds  from  the  north  are 
converging  into  the  western  side  of  the  storm.  The  typhoon  begins  to  accelerate  to  the 
northeast  as  the  steering  flow  around  the  periphery  of  the  STR,  which  appears  to  have 
intensfied  during  the  last  1 2  h.  Again,  analyses  of  intensification  of  this  system  are  sus¬ 
pect  due  to  the  lack  of  observations.  A  significant  aspect  of  this  analysis  is  that  the  re¬ 
curvature  of  Nelson  is  continuing  due  to  the  interaction  with  the  midlatitude  trough, 
even  though  the  two  systems  remain  distinct  features  separated  by  the  col  area. 

6.  1200  UTC  6  October  1988 

Nelson,  located  at  25.3°N,  128. 9°E  with  winds  of  115  kt,  continues  to  weaken, 
move  northeastward  (at  12  kt)  and  accelerate.  The  most  significant  feature  is  that  the 
tropical  cyclone  is  clearly  to  the  north  of  the  ridge  line  connecting  the  two  highs.  This 
indicates  that  Nelson  is  definitely  embedded  in  the  westerlies.  This  feature  is  obvious  in 
the  300  mb  analysis  and  becomes  apparent  later  in  the  500  mb  analysis.  The  upper-level 
col  (Fig.  7a)  between  the  tropical  cyclone  and  the  midlatitude  trough  continues  to  drift 
to  the  east.  The  downstream  winds  to  the  south  of  Japan  still  provide  a  good  outflow, 
and  Nelson  maintains  typhoon  strength  even  as  it  propagates  to  the  north  over  lower 
ocean  temperatures. 

The  500  mb  analysis  has  Nelson  embedded  in  the  westerlies  (Fig.  7b).  Although 
the  trough  north  of  the  storm  seems  to  have  weakened,  the  winds  between  the  trough 
and  Nelson  remain  strong. 
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7.  7-8  October  1988  (Demise  of  Typhoon  Nelson) 

Nelson  continues  to  weaken  as  it  accelerates  northeastward.  The  typhoon  loses 
its  tropical  characteristics  and  transitions  to  an  extratropical  system  southeast  of  Tokyo, 
Japan  at  1500  UTC  8  October  1988.  By  this  time,  the  storm  was  moving  at  a  speed  of 
28  kt.  Although  it  is  becoming  extratropical,  Nelson  maintains  winds  of  typhoon  in¬ 
tensity.  The  final  warning  on  Nelson  was  issued  on  the  1800  UTC  8  October.  Remnants 
of  the  storm  were  visible  on  satellite  imagery  for  two  days  after  the  final  warning. 

Typhoon  Nelson  was  the  only  supertyphoon  of  the  1988  typhoon  season.  For¬ 
tunately,  the  track  to  the  west-northwest  and  subsequent  recurvature  to  the  northeast 
kept  the  storm  away  from  populated  areas  and  major  damage  was  avoided. 

As  mentioned,  it  is  extremely  difficult  to  analyze  the  interaction  of  Typhoon 
Nelson  with  the  extratropical  cyclone  using  conventional  synoptic  data.  Higher  resol¬ 
ution  data  are  required  to  accurately  depict  these  small-scale  features.  NORAPS  ana¬ 
lyses  and  forecasts  will  be  used  to  describe  these  features  because  the  80  km  resolution 
should  be  sufficient  in  determining  the  interaction  during  recurvature  of  Nelson.  Of 
course,  analysis  of  model  data  is  only  reasonable  if  the  model  provides  a  good  forecast 
of  the  storm  evolution.  Data,  such  as  might  be  obtained  during  the  tropical  cyclone 
motion  field  experiment  during  1990  (Abbey  and  Elsberry  1989),  will  be  necessary  to 
prove  that  the  model  description  is  valid. 

B.  FORECAST  OF  TYPHOON  NELSON 

If  the  NORAPS  fields  are  to  be  used  to  infer  the  IPV  interactions  during  recurva¬ 
ture,  it  is  essential  that  the  NORAPS  accurately  predict  the  movement  of  Typhoon 
Nelson.  According  to  Hodur  (19S8),  N'ORAPS  demonstrated  skill  in  predicting  the  path 
of  storms  during  the  1988  typhoon  season,  especially  for  storms  poleward  of  20°N.  For 
those  storms  equatorward  of  20°N,  NORAPS  had  a  rightward  bias.  Hodur  (1988)  states 
the  major  reason  for  tropical  cyclone  forecast  failures  in  NORAPS  during  1988  was  due 
to  the  inability  of  the  model  to  handle  moisture  and  its  interaction  with  radiation. 
However,  Typhoon  Nelson  was  selected  for  this  study  because  it  did  interact  with  a 
midlatitude  trough  and  because  NORAPS  did  a  reasonably  good  job  in  forecasting  the 
recurvature  to  the  northeast. 

The  forecasters  at  JTWC  use  various  tools  in  forecasting  tropical  cyclone  move¬ 
ment:  Climatology  (CLIM),  Half  Persistence  and  Climatology  (HPAC),  two  analog 
programs  (RECR  and  TOTL),  the  One-Way  Tropical  Cyclone  Model  (OTCM),  Ex¬ 
trapolation  (XTRP)  and  Climatology  and  Persistence  (CLIPER).  A  discussion  of  these 
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aids  can  be  found  in  the  ATCR,  1988.  The  various  forecast  tracks  for  Typhoon  Nelson 
are  shown  in  Fig.  8.  Best  track  positions  from  1200  UTC  1  October  until  0000  UTC  4 
October,  and  the  subsequent  forecast  positions  at  24,  48  and  72  h,  are  shown. 

The  majority  of  the  forecast  schemes  accurately  forecast  the  recurvature  of  Nelson. 
The  largest  differences  among  the  schemes  are  related  to  how'  fast  and  how  far  west 
Nelson  would  propagate.  The  JTWC  and  NORAPS  forecasts  have  Nelson  farther  to  the 
east  than  the  other  schemes,  especially  at  72  h.  The  center  of  the  storm  in  the  NORAPS 
forecast  was  assumed  to  be  at  the  center  of  the  maximum  in  the  IPV  fields  presented  in 
the  next  section. 

The  post-storm  analysis  (best  track)  at  1200  UTC  4  October  1988  had  Typhoon 
Nelson  located  at  1S.5°N,  124. 7°E,  whereas  NORAPS  had  the  storm  located  at  20.3°N, 
126. 1°E.  This  200  km  error  to  the  east  indicates  that  the  model  recurved  the  storm  too 
soon.  By  0000  UTC  5  October  198S  Nelson  began  to  recurve  to  the  northeast  and  was 
located  at  21.3°N.  124. 4°E.  The  NORAPS  forecast  position  was  at  21.3°N,  126.4°E. 
Although  the  storm  actually  remained  relatively  stationary  in  the  last  12  h,  the  NORAPS 
forecast  has  the  storm  far  to  the  east  because  of  the  early  recurvature.  Nelson  has  begun 
its  propagation  to  the  northeast  as  indicated  by  the  1200  UTC  5  October  19S8  position 
at  22.5°N,  125.1°E.  Although  NORAPS  accurately  predicted  a  track  to  the  northeast, 
the  forecast  position  of  23.1°N,  126.5°E  still  shows  an  eastward  bias.  Nelson  continues 
to  move  to  the  northeast  and  at  0000  UTC  6  October  1988  it  is  located  at  24.0°N, 
126.7°E.  The  48-h  NORAPS  forecast  has  Nelson  positioned  at  24.5°N,  127. 7°E.  Al¬ 
though  there  is  still  an  eastern  bias  in  the  model  forecast,  the  position  error  is  only  148 
km.  At  1200  UTC  6  October  1988  Nelson  is  located  at  25.3°N,  128. 9°E  and  the 
NORAPS  position  is  at  26.8°N,  128.9°E.  Although  the  easterly  bias  in  the  previous 
forecast  periods  has  been  eliminated,  the  model  has  continued  to  accelerate  the  storm 
farther  to  the  north.  For  the  final  forecast  period  {ending  at  0000  UTC  7  October  1988), 
Nelson  is  located  at  26.6°N,  131.1°E,  whereas  the  NORAPS  position  is  at  29.6°N, 
133.6°E.  The  model  has  pushed  the  storm  too  far  to  the  northeast  with  a  72-h  error  of 
388  km. 

Overall,  NORAPS  did  an  excellent  job  in  forecasting  the  movement  of  Typhoon 
Nelson.  As  mentioned  above,  the  relatively  stationary  position  of  the  large-scale  fea¬ 
tures  that  control  the  typhoon  movement  probably  contributes  to  the  accurate  forecast. 
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Figure  8.  Tracks  of  Typhoon  Nelson  forecast  by  various  aids  where  XTRP(X), 
TOTL(T),  HPAC(H),  CLIPPER(C),  RECR(R),  OTCM(O),  JTWC(J),  CLIM(M), 
NORAPS(N)  and  Best  Track(B).  The  letters  along  each  track  indicate  the  positions 
at  24,  48  and  72  h  after  the  initial  time  of  0000  UTC  4  Oct  I98S.  Earlier  storm  po¬ 
sitions  with  a  12-h  interval  are  indicated  by  +. 
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IV.  DESCRIPTION  OF  TYPHOON  NELSON  IN  TERMS  OF  IPV. 


The  synoptic  analysis  of  Typhoon  Nelson  using  NOGAPS  300  and  500  mb  height 
fields  indicates  that  the  interaction  between  the  midlatitude  and  tropical  cyclones  occurs 
on  horizontal  scales  that  are  too  small  to  be  accurately  represented  by  synoptic-scale 
data.  The  major  problem  is  the  deficiency  in  horizontal  and  vertical  resolution  due  to 
the  lack  of  observations  in  the  vicinity  of  the  tropical  cyclone.  By  using  model  fields 
with  80  km  horizontal  resolution  and  21  level  vertical  resolution,  these  "gaps"  in  data 
are  filled  with  dynamically  consistent  model  data.  The  analysis  of  isentropic  potential 
vorticity  (IPV)  fields  in  Typhoon  Nelson  illustrates  two  interesting  features.  First,  the 
vertical  and  horizontal  structure  of  the  bogus  vortex  is  analyzed  to  determine  the  effect 
it  may  have  on  the  NORAPS  forecast.  Second,  the  IPV  advection  pattern  associated 
with  the  storm  center  and  other  synoptic  features  are  examined  in  relation  to  storm 
propagation. 

A.  HORIZONTAL  AND  VERTICAL  STRUCTURE  OF  THE  BOGUS  STORM 

Hodur  (1988)  states  that  the  bogus  vortex  used  in  NORAPS  is  much  larger  than 
typical  tropical  cyclones  and  this  may  degrade  the  track  forecast,  especially  in  low  lati¬ 
tude  forecasts.  The  large  radial  extent  of  the  bogus  storm  is  seen  in  the  horizontal  rep¬ 
resentation  of  Typhoon  Nelson  in  terms  of  IPV  (Fig.  9a).  The  center  of  Typhoon 
Nelson  at  310  K  corresponds  to  the  maximum  value  of  IPV  because  the  relative  vorticity 
is  a  maximum  inside  the  radius  of  maximum  winds.  Large  values  of  IPV  encircle  the 
storm  center,  with  highest  values  to  the  extreme  east  and  west  of  the  storm.  The  radial 
width  of  the  IPV  associated  with  the  storm  is  surprisingly  broad  for  the  low  levels  of  a 
tropical  cyclone.  Moreover,  the  horizontal  scale  of  the  cyclonic  winds  also  appears  to 
be  unrealistically  large  as  the  bogus  vortex  dominates  the  western  North  Pacific  area. 

The  bogus  storm  is  clearly  depicted  on  the  320  K  surface  (Fig.  9b).  Although 
isolines  of  IPV  still  encircle  the  storm  center,  the  vortex  is  not  circularly  symmetric  at 
larger  radii  because  a  larger  node  of  IPV  extends  southward  and  a  maximum  of 
3.5  x  106mKs'3mb*'  is  found  just  north  of  Taiwan.  Another  secondary  maximum  is 
found  in  the  southwest  corner.  The  large  gap  between  the  tropical  storm  and  the  max¬ 
imum  IPV  associated  with  midlatitude  systems  also  is  evident. 

The  upper-tropospheric  IPV  representation  of  the  initial  bogus  storm  is  shown  in 
Fig.  10a.  Maximum  IPV  values  are  found  at  the  340  K  level.  A  strong  westerly  jet  is 


21 


Figure  9.  IPV  at  0000  UTC  4  Oct  1988  on  (A)  310  K  and  (B)  320  K  surfaces. 
Wind  speeds  are  proportional  to  the  length  of  the  vectors  and  IPV  contours  are 
1  x  lO^rnKs^nib-’  .  Heavy  dashed  line  is  the  800  km  radius. 
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depicted  along  the  northern  boundary  with  a  IPV  dipole  of  high  (low)  values  on  the 
cyclonic  (anticyclonic)  side  to  the  northeast  of  the  storm  center.  A  ring  of  low  IPV 
values  encircles  the  storm  and  the  values  are  actually  negative  to  the  west  and  southwest 
of  the  vortex.  These  minima  are  on  the  boundary  between  cyclonic  flow  in  the  storm 
and  the  anticyclonic  outflow.  An  IPV  maximum  is  found  northwest  of  the  vortex  where 
the  anticyclonic  flow  on  the  edge  of  the  vortex  interacts  with  the  westerly  wind  maxi¬ 
mum  to  the  north.  This  feature  is  also  found  higher  in  the  troposphere  on  the  350  K 
surface  (Fig.  10b).  In  the  350  K  surface,  another  maximum  of  IPV  larger  than  the  bo¬ 
gus  vortex  is  located  just  west-northwest  of  the  vortex  center.  The  values  of  IPV  near 
the  center  of  the  vortex  are  smaller  then  at  the  340  K  surface  as  the  cyclonic  flow  around 
the  center  is  not  as  strong.  However,  large  negative  values  of  IPV  are  found  to  the  south 
of  of  the  vortex  center  along  the  edge  of  a  band  of  strong  easterly  flow.  In  summary, 
this  upper-level  IPV  representation  of  the  tropical  cyclone  and  the  midlatitude  flow 
would  seem  to  indicate  a  strong  interaction  is  already  occurring  with  southerly  and 
southwesterly  flow  to  the  northwest  of  Nelson. 

The  vertical  structure  of  these  features  relative  to  the  vertical  structure  of  the  bogus 
vortex  is  indicated  by  a  vertical  cross-section  from  29°\,  119°E  to  12°N,  132°E  in 
Fig.  11a.  The  lower  boundary  corresponds  to  the  300  K  surface  and  the  upper  bound¬ 
ary  is  along  the  350  K  surface.  Limiting  the  upper  boundary  to  350  K  removes  the  ex¬ 
tremely  large  IPV  values  associated  with  the  stratospheric  stratification.  These 
stratospheric  features  are  not  considered  relevant  for  this  study.  Cyclonic  flow  around 
the  vortex  is  represented  as  positive  negative  wind  values  to  the  right  left  of  the  storm 
center,  which  is  represented  by  the  zero  isoline.  The  warm  core  structure  of  the  bogus 
vortex  is  indicated  by  the  downward  depression  of  the  350  K  surface  over  the  vortex 
center. 

In  Fig.  I  la  the  upper-tropospheric  IPV  maximum  described  above  is  located  to  the 
left  (northwest)  of  the  vortex.  This  feature  extends  into  the  lower  troposphere  with  a 
secondary  maximum  at  approximately  750  mb.  A  narrow  region  of  negative  IPV  values 
extends  downward  to  about  400  mb  between  the  IPV  maximum  to  the  northwest  and 
the  maximum  associated  with  the  tropical  cyclone.  This  negative  IPV  region  is  on  the 
anticyclonic  shear  side  of  the  zero  isoline  in  the  upper  tropospheric  outflow  region. 
Another  region  of  negative  values  of  IPV  located  southwest  of  the  vortex  is  illustrated 
by  the  cross-section  in  Fig.  lib.  An  interesting  feature  in  both  cross-sections  is  the 
vertical  coupling  of  the  upper  and  lower  vortex  structures.  In  both  displays,  the  zero 
velocity  isoline  shows  very  little  horizontal  displacement.  The  large  radial  size  of  the 
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bogus  vortex  is  also  evident  in  the  cross-sections  as  the  positive  IPV  values  around  the 
vortex  are  extremely  broad.  This  is  another  indication  that  the  bogus  vortex  is  too  large 
(Hodur  19S8,  1989). 

After  a  12  h  integration,  the  large  radial  extent  of  the  bogus  vortex  is  still  apparent 
(Fig.  12a).  The  bogus  vortex  has  become  more  asymmetric  as  a  lobe  of  low  IPV  has 
developed  just  southwest  of  the  center.  The  cyclonic  circulation  is  still  extensive  with  a 
good  inflow  from  the  southwest.  At  320  K  (Fig.  12b),  the  vortex  has  become  more 
compact.  However,  the  circulation  around  the  storm  also  includes  asymmetries.  These 
asymmetries  may  be  due  to  adjustment  in  the  model  to  compensate  for  the  large  warm 
core  system  represented  by  the  bogus  vortex.  In  the  upper  troposphere  (Fig.  1 3a),  the 
IPV  field  around  the  vortex  has  become  extremely  compact  and  oval  in  shape.  The  IPV 
field  is  oriented  to  north-south  and  with  a  significant  wavenumber  2  component.  This 
is  a  change  from  the  initial  conditions  when  the  wind  field  was  more  circular.  Negative 
values  of  IPV  are  located  to  the  southwest  and  southeast  of  the  storm  center.  An  area 
with  a  large  IPV  maximum  is  found  to  the  northeast  of  the  storm  center.  These  values 
are  also  evident  in  Fig.  13b.  At  this  level,  the  cyclonic  flow  around  the  vortex  has  di¬ 
minished. 

An  interesting  feature  that  occurs  after  12  h  is  the  vertical  displacement  of  the 
maximum  value  of  IPV  associated  with  the  vortex  center.  In  the  initial  bogus,  the 
lower-level  and  upper-level  features  were  well  coupled.  After  a  12  h  integration,  the  IPV 
maxima  are  not  vertically  stacked  (compare  Fig.  12b  to  Fig.  13a).  The  displacement  to 
the  south  of  the  maximum  value  at  upper  levels  is  readily  apparent  in  Fig.  14a.  A 
"wobble"  has  been  introduced  into  the  vortex  structure.  This  uncoupling  of  the  upper 
and  lower  level  may  be  due  to  the  inability  of  the  model  to  sustain  the  bogus  vortex. 
The  presence  of  negative  IPV  on  the  periphery  of  the  bogus  vortex  may  have  led  to 
inertial  instability  (Ciesielski  et  al.  1989),  which  displaces  high  values  of  PV  along  the 
periphery  of  the  vortex.  This  unstable  "wobble"  is  also  seen  in  Fig.  14b  as  the  low-level 
PV  maximum  is  to  the  southeast  of  the  upper-level  maximum.  Another  interesting  fea¬ 
ture  in  both  cross-sections  is  the  difference  of  the  maximum  winds  on  the  eastern  side 
of  the  vortex  as  compared  to  those  on  the  western  side.  Since  this  difference  in  velocity 
is  due  to  the  superposition  of  a  vortex  on  the  advecting  flow,  it  should  indicate  propa¬ 
gation  to  the  north. 

Clear  indications  of  stability  in  the  bogus  vortex  occur  after  36  h  (Fig.  15a).  To 
accommodate  the  movement  of  the  storm,  the  grid  has  been  displaced  to  the  north.  The 
most  noticeable  feature  is  the  reduction  in  the  size  of  the  vortex  compared  to  earlier 
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Figure  11.  Cross-section  on  00  UTC  4  Oct  19S8  from  (A)  northwest  to  southeast 
and  (B)  southwest  to  northeast.  Potential  vorticity  { lO'mKs  hnb  dashed)  and 
tangential  wind  components  (m.s,  solid)  arc  shown  between  300  K  (lower  boundary) 
and  350  K  (upper  boundary).  Orientation  of  cross-sections  shown  in  Fig.  10. 
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Figure  14.  As  in  Fig.  11,  except  for  12  UTC  4  Oct  1988  for  (A)  northwest  to  the 
southeast  and  (B)  southwest  to  the  northeast. 


times.  The  IPV  maximum  near  the  vortex  center  is  more  compact,  and  lobes  of  high 
IPV  around  the  periphery’  of  the  vortex  arc  not  as  numerous.  However,  the  maximum 


value  of  PV  associated  with  the  vortex  center  has  almost  doubled  in  36  h.  The  model 
evidently  has  adjusted  to  the  dynamical  processes  associated  with  the  tropical  cyclone 
and  the  1PV  maximum  now  depicts  the  intense  diabatic  heating  found  near  the  storm 
center  in  the  model. 

High  values  of  IPV  along  the  periphery  of  the  vortex  are  still  apparent  after  36  h  in 
Fig.  15b,  which  may  indicate  that  the  model  is  still  trying  to  adjust  to  the  initial  wind 
and  mass  fields.  Although  the  IPV  maximum  at  this  level  has  also  increased,  the  vortex 
is  not  as  symmetric.  The  lobes  of  high  IPV  to  the  east  of  the  vortex  may  be  associated 
with  the  greater  diabatic  heating  due  to  the  advection  of  warm  moist  air  from  the 
southwest  quadrant. 

In  the  upper  troposphere,  the  vortex  circulation  after  36  h  is  very  compact 
(Fig.  16a).  Small  negative  values  of  IPV  are  found  in  the  anticyclonic  cell  to  the  east 
of  the  center.  This  feature  is  further  enhanced  at  350  K  (Fig.  16b)  as  negative  values 
of  IPV  are  found  to  the  east  and  west  of  the  vortex.  The  lobe  of  cyclonic  IPV  to  the 
northwest  of  the  vortex  represents  an  approaching  midlatitude  wave.  The  vertical  extent 
of  this  feature  is  evident  in  Fig.  17a.  The  maximum  of  IPV  associated  with  this  short¬ 
wave  extends  downward  into  the  midtroposphere  to  approximately  4S0  mb.  This 
cross-section  also  illustrates  that  the  vortex  has  become  much  more  compact  after  36 
h.  The  lower-level  and  upper-level  maxima  of  IPV  are  much  better  coupled.  The  stable 
structure  of  the  vortex  is  also  apparent  in  Fig.  17b,  which  again  has  very  little  vertical 
tilt.  Since  the  "wobble"  in  the  vortex  has  diminished,  the  model  has  adjusted  the  mass 
and  wind  fields  to  the  diabatic  heating  associated  with  the  vortex.  If  this  structure  better 
represents  the  "model  dynamics"  of  the  tropical  cyclone,  the  bogus  vortex  that  is  inserted 
into  the  initial  fields  should  resemble  this  structure  rather  than  the  broad  vortex  at  0000 
UTC  4  October. 

B.  ADVECTION  OF  IPV  FIELDS 

Forecasters  at  the  Joint  Typhoon  Warning  Center  (JTWC)  make  use  of  various 
objective  aids  to  guide  them  in  the  forecast  of  tropical  cyclones.  Since  only  the  track 
forecast  and  not  the  predicted  fields  are  transmitted  to  JTWC,  the  forecaster  has  little 
opportunity  to  assess  the  dynamics  involved  in  the  model  forecast.  Without  dynamical 
reasoning,  the  likely  accuracy  of  the  model  forecast  is  difficult  to  assess.  If  the  JTWC 
forecaster  knew  what  dynamical  features  contributed  to  the  track  forecast,  he,  she  would 
be  able  to  assess  the  synoptic  realism  and  monitor  the  relevant  synoptic  features  during 
the  course  of  the  forecast. 
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Figure  15.  As  in  Fig.  9,  except  for  12  UTC  5  Oct  1988  at  (A)  310  K  and  (B)  320 

K. 
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Figure  17.  As  in  Fig.  11,  except  for  12  UTC  5  Oct  for  (A)  northwest  to  southeast 
and  (B)  southwest  to  northeast. 
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The  purpose  of  this  section  is  to  illustrate  the  usefulness  of  a  IPV  advection  field 
representation,  and  specifically  how  this  may  aid  a  forecaster  in  determining  the  con¬ 
tributions  to  storm  movement  due  to  self-advection,  advection  due  to  basic  flow  and  the 
advection  resulting  from  adjacent  synoptic  features.  Two  theta  surfaces  (310  K  and  320 
K)  are  selected  to  represent  the  mid-tropospheric  steering  flow.  The  two  upper-level 
surfaces  (340  K  and  350  K)  are  presented  for  the  initial  period  and  after  12  h  to  show 
the  effect  of  the  bogus  vortex  on  the  IPV  advection  fields.  The  IPV  advection  fields  are 
just  one  of  the  output  fields  generated  by  NORAPS  that  may  be  useful  to  the  JTWC 
forecaster. 

At  0000  UTC  4  October  1988,  a  large  region  of  positive  IPV  advection  (IPVA)  is 
located  to  the  northwest  of  the  storm  center  (Fig.  18a).  The  positive/negative  IPVA 
dipole  associated  with  the  flow  around  the  vortex  is  the  dominating  feature  at  this  time. 
A  vector  from  the  negative  IPVA  center  to  the  positive  IPVA  center  indicates  the  in¬ 
stantaneous  IPVA  tendencies,  and  suggests  a  northwestward  track.  Since  no  other  sig¬ 
nificant  IPVA  centers  are  found  near  the  vortex,  it  is  expected  that  this  motion  tendency 
will  continue  for  some  time.  A  similar  lobe  configuration  is  found  at  320  K  (Fig.  18b), 
except  both  lobes  are  now  asymmetric  with  a  northeast-southwest  elongation.  However, 
the  orientation  of  the  IPVA  centers  with  respect  to  the  vortex  still  indicates  movement 
to  the  northwest.  As  in  the  preceding  figure,  the  advection  pattern  associated  with  the 
vortex  is  the  dominant  feature. 

Maxima  of  IPV  advection  that  are  not  associated  with  the  storm  center  are  found 
in  the  upper  troposphere  (Fig.  1 9a).  Lobes  of  positive  and  negative  IPVA  are  located 
northwest  of  the  storm  center.  Although  these  lobes  might  appear  to  represent  a  mid¬ 
latitude  wave,  they  are  actually  a  result  of  the  extremely  broad  bogus  vortex.  The 
advection  pattern  near  the  center  is  asymmetric  as  both  lobes  are  elongated  to  the 
southwest.  Connecting  the  maximum  values  of  negative  and  positive  IPVA  again  im¬ 
plies  movement  to  the  northwest.  The  effect  of  the  over-sized  bogus  is  also  apparent  in 
Fig.  19b.  Maxima  and  minima  values  of  IPVA  are  located  along  the  outer  periphery 
of  the  bogus  vortex.  Thus,  the  bogus  vortex  distorts  the  IPVA  fields,  which  may  lead 
to  an  inaccurate  representation  of  the  motion  processes  at  this  level. 

After  a  12  h  integration,  little  change  occurs  in  the  lower  tropospheric  advection 
pattern  (Fig.  20a).  The  positive  (negative)  IPVA  center  is  still  to  the  northwest 
(southeast)  of  the  vortex  center.  The  IPV  advection  pattern  associated  with  the  vortex 
remains  the  dominant  feature.  At  320  K  (Fig.  20b),  the  advection  patterns  near  the 
vortex  have  become  more  distorted.  The  negative  lobe  is  extended  farther  to  the  south 
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Figure  18.  IPV  advection  at  00  UTC  4  Oct  1988  Tor  (A)  310  K  and  (B)  320  K. 
Positive/negative  IPV  advection  (contour  intervals  0.7  for  (A)  and  0.6  for  (B) 
KPKs^mb-1)  are  indicated  by  II, 'L,  the  position  of  Typhoon  Nelson  is  indicated  by 
the  dot  and  the  wind  speeds  are  propotional  to  the  vector  lengths. 
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Figure  19.  As  in  Fig.  18,  except  for  00  U'l'C  4  Oct  1988  at  (A)  340  K  and  (B)  350 
K.  Contour  intervals  are  0.7  for  (A)  and  4  for  (B) 
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as  two  minima  occur.  Nevertheless,  self-advection  is  likely  the  key  process  influencing 
the  propagation  of  the  storm. 

At  340  K  (Fig.  21a),  other  1PV  advection  features  are  apparent.  Significant  positive 
and  negative  IPVA  centers  are  located  to  the  northeast  of  the  vortex.  Later  analyses 
(not  shown)  indicate  that  these  lobes  are  transient  in  nature  and  therefore  are  of  little 
significance  in  storm  propagation..  The  two  lobes  near  the  storm  center  have  become 
extremely  distorted  as  each  lobe  now  has  two  centers.  It  would  be  difficult  to  determine 
the  vortex  motion  from  this  IPVA  pattern.  The  "wobble”  of  the  upper-level  maxima 
described  earlier  is  probably  related  to  this  complex  IPVA  pattern. 

The  effects  of  the  bogus  vortex  on  the  IPVA  pattern  at  350  K  has  decreased  in  12 
h  (Fig.  21b).  Although  the  pattern  still  shows  some  distortion,  the  lobes  of  positive  and 
negative  IPVA  are  not  as  numerous  as  is  in  the  initial  pattern  (Fig.  19b).  The  dominant 
upper-level  IPVA  near  the  vortex  shows  very  little  change  from  the  preceding  analysis. 

Significant  advection  maxima  other  than  those  associated  with  the  vortex  are  first 
detected  at  1200  L’TC  5  October  I9S8  (Fig.  22a).  A  dipole  of  positive  and  negative 
IPVA  centers  is  found  to  the  northeast  of  the  vortex  and  is  probably  due  to  the  large- 
scale  synoptic  flow.  The  dipole  pattern  associated  with  the  vortex  indicates  a  north- 
northwest  propagation  at  this  time,  which  agrees  with  the  predicted  track.  At  320  K 
(Fig.  22b).  the  IPVA  maxima  are  more  distorted.  The  maximum  value  of  IPVA  is  al¬ 
most  directly  north  of  the  storm  center  and  two  minima  occur,  one  south  of  the  center 
and  the  other  to  the  southeast.  Consequently,  the  orientation  of  the  propagation  vector 
is  not  well  established  at  this  level.  However,  the  extension  of  the  positive  IPVA  toward 
the  north  and  northeast  suggests  continued  recurvature.  In  the  terminology  of  Sherman 
(1988),  this  would  be  a  "primary  node"  toward  which  the  storm  would  propagate. 

By  0000  LTC  6  October  1988,  the  advection  field  has  become  more  complicated 
(Fig.  23a).  The  centers  of  positive  and  negative  IPVA  near  the  vortex  are  still  apparent. 
However,  other  centers  have  developed  to  the  north  and  northwest  of  the  vortex.  The 
dipole  to  the  northwest  is  due  to  a  propagating  shortwave.  The  negative  lobe  to  the 
north  contributes  to  a  "wavetrain"  pattern  oriented  along  the  general  synoptic  flow. 
This  maximumminimum  pattern  to  the  northeast  may  be  a  key  indication  of  the  recur¬ 
vature  of  the  storm.  This  pattern  is  also  apparent  at  320  K  (Fig.  23b).  The 
positive,  negative  wavetrain  of  IPVA  to  the  northeast  of  the  vortex  may  be  a  precursor 
of  the  subsequent  storm  movement  to  the  northeast.  Since  this  is  only  a  single  case 
study,  the  general  applicability  of  this  feature  must  be  determined  from  many  more 
cases. 
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Figure  20.  As  in  Fig.  18,  except  for  12  UTC  4  Oct  1988  at  (A)  310  K  and  (B)  320 
K.  Contour  intervals  are  0.6  for  (A)  and  0.6  for  (B). 
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Figure  21.  As  in  Fig.  19,  except  at  12  UTC  4  Oct  1988  for  (A)  340  K  and  (B)  350 
K.  Contour  intervals  are  1  for  (A)  and  5  lor  (B). 


Figure  22.  As  in  Fig.  18,  except  at  12  UTC  5  Oct  1988  for  (A)  310  K  and  (B)  320 
K.  Contour  intervals  are  1  for  (A)  and  1  for  (B). 


40 


Figure  23.  As  in  Fig.  18,  except  at  00  UTC  6  Oct  1988  for  (A)  310  K  and  (B)  320 
K.  Contour  intervals  are  1  for  (A)  and  2  for  (B). 


At  1200  UTC  6  October  1988,  the  centers  of  positive  and  negative  IPVA  are  still 
found  to  the  northeast  of  the  storm  (Fig.  24a)  as  in  the  pattern  24  h  previously.  The 
positive  and  negative  centers  of  IPVA  associated  with  the  short  wave  are  also  readily 
apparent  in  this  figure.  This  wave  continues  to  propagate  to  the  southeast.  At  320  K 
(Fig.  24b),  the  short-wave  and  the  large-scale  feature  to  the  northeast  of  the  vortex  are 
not  well  represented. 

At  the  final  time  (Fig.  25a),  the  configuration  of  the  310  K  pattern  with  repect  to 
the  vortex  has  not  changed  significantly  in  12  h.  A  large  center  of  negative  IPVA  is 
situated  to  the  northeast  of  the  storm  center.  Although  the  tropical  cyclone  is  deeply 
embedded  in  the  large-scale  flow,  the  dipole  associated  with  the  vortex  shows  very  little 
distortion.  The  shortwave  to  the  west  of  the  storm  is  not  well  organized  as  it  continues 
to  propagate  to  the  southeast.  A  similar  orientation  of  the  IPVA  is  found  at  320  K 
(Fig.  25b).  The  large-scale  orientation  to  the  northeast  is  more  apparent  as  the  negative 
IPVA  center  interacts  with  the  IPVA  dipole  associated  with  the  tropical  cyclone. 

In  summary,  significant  changes  are  revealed  in  the  IPVA  fields  as  the  storm  evolves 
from  the  initial  analysis  to  the  forecast  fields  after  72  h.  Self-advection  was  initially  the 
major  influence  on  storm  movement,  and  adjacent  advective  features  only  appear  at  later 
times.  The  significance  of  these  features  and  their  role  in  storm  propagation  can  be  de¬ 
termined  only  by  further  study. 

The  transmission  of  these  IPV  fields  along  with  the  NORAPS  track  forecast  may 
contribute  to  the  understanding  of  the  dynamical  processes  involved  in  storm  propa¬ 
gation.  Such  fields  could  contribute  to  the  decision  of  whether  the  recurvature  forecast 
of  the  tropical  cyclone  is  synoptically  consistent  and  thus  should  be  accepted. 
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Figure  24.  As  in  Fig.  18.  except  at  12  UTC  6  Oct  1988  for  (A)  310  K  and  (B)  320 
K.  Contour  intervals  are  1  for  (A)  and  2  for  (B). 


Figure  25.  As  in  Fig.  18,  except  at  00  U  1C  7  Oct  1988  for  (A)  3)0  K  and  (B)  320 
K.  Contour  intervals  are  2  for  (A)  and  2  for  (B). 
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V.  CONCLUSIONS 


The  purpose  of  this  study  is  to  determine  the  feasibility  of  using  model-generated 
isentropic  potential  vorticity  (IPV)  fields  to  aid  in  the  forecast  of  tropical  cyclone 
movement.  The  IPV  fields  are  from  the  NORAPS  forecast  of  Typhoon  Nelso.,  m  the 
period  0000  UTC  4  October  to  0000  UTC  7  October  1988.  The  fields  are  analyzed  every 
12  h. 

Various  researchers  have  made  use  of  the  distribution  and  structure  of  IPV  fields  in 
their  attempt  to  explain  midlatitude  phenomena  such  as  cyclogenesis  and  explosive 
cyclogenesis.  The  major  deficiency  in  these  studies  is  the  lack  of  horizontal  and  vertical 
resolution  due  to  the  sparsity  of  observations  in  the  area  of  interest.  This  is  also  a  major 
concern  in  tropical  analysis  where  few  rawinsondes  exist.  To  overcome  this  obstacle, 
researchers  have  used  numerical  models.  However,  most  operational  numerical  models 
do  not  have  adequate  vertical  resolution  to  describe  key  IPV  features.  NORAPS  was 
selected  for  this  study  because  its  vertical  (21  levels)  and  horizontal  (SO  km)  resolution 
should  be  sufficient  to  adequately  describe  these  features. 

A  synoptic  analysis  of  Typhoon  Nelson  using  NOGAPS  300  and  500  mb  height 
fields  is  used  first  to  determine  if  Nelson  interacted  with  a  midlatitude  cyclone.  However, 
the  interaction  appears  to  occur  on  horizontal  scales  too  small  to  be  accurately  repres¬ 
ented  by  synoptic-scale  data.  The  motivation  for  using  model-generated  fields  is  that 
these  "gaps"  in  data  are  avoided. 

The  analysis  of  IPV  fields  in  Typhoon  Nelson  illustrates  how  the  over-sized  bogus 
vortex  distorts  the  IPV  fields  and  possibly  degrades  the  storm  forecast.  The  model  evi¬ 
dently  responds  to  negative  IPV  values  on  the  periphery  of  the  large  bogus  via  inertial 
instability.  An  uncoupling  of  the  upper  and  lower  levels  may  be  due  to  the  inability  of 
the  model  to  sustain  the  bogus  vortex.  The  over-sized  bogus  may  lead  to  a  degradation 
in  the  model  forecast  as  it  takes  the  model  approximately  36  h  to  adjust  to  the  initial 
imbalance  in  the  mass  and  wind  fields. 

The  IPV  advection  patterns  are  examined  to  determine  if  the  calculations  of  these 
fields  may  aid  the  forecaster  in  isolating  the  key  features  governing  tropical  cyclone 
movement.  Significant  changes  in  the  1PVA  fields  are  found  as  the  storm  evolved  from 
the  initial  analysis  to  the  forecast  field  72  h  later.  Initially,  self-advection  due  to  the 
orientation  of  the  positive  and  negative  dipoles  near  the  storm  center  is  the  dominating 
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factor  in  the  storm  movement.  After  a  36  h  integration,  other  adjacent  advective  fea¬ 
tures  become  apparent. 

In  future  research,  IPV  fields  for  additional  tropical  cyclones  undergoing  recurvature 
should  be  analyzed.  The  general  applicability  of  the  interaction  of  IPV  fields  and  their 
apparent  effects  on  storm  propagation  in  this  case  study  then  can  be  determined.  An¬ 
other  illustration  of  the  significance  of  these  interactions  is  to  modify  the  initial 
NORAPS  fields  so  that  the  tropical  cyclone  is  displayed  southward  from  the  initial  po¬ 
sition.  An  examination  of  the  differences  in  the  IPV  interactions  when  the  modified 
storm  does  not  recurve  should  indicate  the  IPV  effects  to  distinguish  between  recurva¬ 
ture  and  non-recurvature. 

An  effort  must  be  made  to  improve  the  bogus  structure.  A  bogus  vortex  that  better 
represents  the  actual  size  of  the  tropical  cyclone  is  needed,  especially  in  the  upper  levels. 
A  better  representation  of  the  vortex  should  reduce  the  time  required  for  the  model  to 
adjust  to  the  imbalance  in  the  mass  and  wind  fields.  A  realistic  bogus  thus  may  improve 
the  early  portion  of  the  model  forecast. 

Finally,  an  effort  must  be  made  to  transmit  dynamical  fields,  such  as  IPV,  to  the 
forecasters.  An  enhancement  of  the  forecaster's  dynamical  knowledge  will  lead  to  im¬ 
provements  in  tropical  cyclone  forecasts. 
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APPENDIX  A.  NAVY  OPERATIONAL  REGIONAL  ATMOSPHERIC 

PREDICTION  SYSTEM 

NORAPS  was  developed  by  Dr.  Rich  Hodur  of  the  Naval  Environmental  Prediction 
Research  Facility  (NEPRF).  The  purpose  of  NORAPS  is  to  provide  a  high  resolution 
forecast  of  up  to  42  h  for  those  areas  which  are  of  importance  to  military  operations. 
Since  1987  NORAPS  has  been  used  as  a  tropical  cyclone  track  prediction  model  in  the 
western  North  Pacific  Ocean.  This  has  been  accomplished  by  artificially  incorporating 
tropical  cyclone(s)  into  the  NORAPS  initial  fields  and  tracking  their  positions  during  the 
forecast  (Hodur  1988). 

The  major  advantage  that  NORAPS  has  over  global  or  hemispheric  models  is  the 
capability  to  resolve  very  small  spatial  features.  Another  key  attribute  of  NO  RAPS  is 
its  flexibility.  NORAPS  is  designed  to  be  used  anywhere  in  the  world  with  a  globally 
relocatable  grid.  The  user  has  the  ability  to  specify  the  dimensions  and  the 
horizontal  vertical  resolution.  An  additional  feature  of  NORAPS  is  the  "terrain  envel¬ 
oping"  concept  in  which  topography  is  calculated  at  a  high  horizontal  resolution  so  as 
to  incorporate  the  effects  of  the  sub-grid  scale  features  into  the  topography  field.  A 
more  indepth  discussion  of  NORAPS  is  provided  by  Hodur  (1982,  1987). 

NORAPS  consists  of  four  major  components:  analysis,  initialization,  forecast  and 
output.  In  the  analysis,  NORAPS  includes  a  regional  analysis  of  wind,  temperature  and 
D  values.  In  the  operational  version  described  by  Hodur  (1987),  first-guess  values  are 
obtained  by  interpolating  Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  fields  with  a  bicubic  spline  to  the  NORAPS  forecast  domain.  Sea-level 
pressure  and  sea-surface  temperatures  obtained  from  the  Fleet  Numerical  Oceanography 
Center  (FN'OC)  are  also  incorporated  into  the  analysis.  These  two  analyses  are  run  in¬ 
dependently  of  NORAPS  and  NOGAPS. 

NORAPS  now  uses  a  regional  multivariate  optimum  interpolation  (01)  scheme  us¬ 
ing  15  pressure  levels  (R.  Hodur  1988).  This  method  is  based  on  the  work  of  Lorenc 
(19S1)  and  includes  a  second-order  autoregression  function  as  the  forecast  error  corre¬ 
lation  function.  The  incorporation  of  the  01  scheme  is  a  departure  from  the  operational 
version  of  NORAPS  as  described  by  Hodur  (1987).  Another  modification  to  the  oper¬ 
ational  version  is  the  use  of  a  new  initialization  phase.  The  research  version  used  in  this 
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study  incorporates  a  nonlinear  vertical  mode  initialization  as  described  by  Bourke  and 
McGregor  (1983). 

Rawinsonde,  aircraft,  ship  and  pibal  data,  as  well  as  satellite  soundings  are  used  in 
the  analysis  procedure.  Hodur  (1987)  states  that  satellite  data  are  used  exclusively  over 
the  ocean.  All  observation  types  are  weighted  equally  in  the  wind  analysis,  whereas 
rawinsonde  data  are  weighted  three  times  larger  than  all  other  data  types  in  the  mass 
analysis. 

The  forecast  component  is  the  main  driving  force  of  NORAPS  and  requires  the 
majority  of  the  computer  time.  NORAPS  uses  the  flux  form  of  the  primitive  equations 
on  a  staggered  grid  scheme  C  as  described  by  Arakawa  and  Lamb  (1977).  The  vertical 
structure  of  the  grid  is  shown  in  Fig.  26  and  the  horizontal  stucture  is  shown  in 
Fig.  27.  All  of  the  prognostic  variables  (u,  v,  q  and  T)  are  located  at  the  middle  of  each 
layer,  and  the  diagnostic  vertical  velocity  is  at  the  sigma  level.  The  model  uses  sigma 
as  the  vertical  coordinate,  which  orients  all  coordinate  surfaces  parallel  to  the  terrain 
surface.  Sigma  is  defined  as  the  pressure  at  any  height  divided  by  the  surface  pressure. 
The  research  version  of  the  forecast  model  uses  21  sigma  levels  (Table  1)  instead  of  the 
12  levels  described  by  Hodur  (1987). 

The  time  integration  scheme  in  the  forecast  model  is  a  split-explicit  technique  de¬ 
veloped  by  Madala  (1981).  Even  though  this  scheme  allows  for  longer  time  steps  for  the 
slower  meteorological  modes,  it  still  has  the  ability  to  predict  all  of  the  gravity  modes. 
The  specific  size  of  the  time  step  is  governed  by  the  computational  stability  criterion  for 
the  horizontal  resolution  selected.  All  of  the  equations  in  the  model  are  solved  using  a 
centered  time  (leapfrog)  and  space  differencing  scheme.  Fourth-order  advection  is  used 
in  all  of  the  equations. 

Hodur  (1987)  states  that  fourth-order  diffusion  is  used  in  all  the  prognostic- 
equations,  except  that  second-order  diffusion  is  used  on  the  first  interior  row.  The  dif¬ 
fusion  is  done  on  the  sigma  levels,  which  has  a  tendency  to  lead  to  spurious  sources  and 
sinks  of  energy  when  applied  to  the  mass  and  moisture  fields.  This  problem  is  alleviated 
by  applying  diffusion  to  deviations  from  the  standard  atmosphere  for  the  temperature, 
moisture  and  pressure.  The  diffusion  coefficient  is  then  computed  in  the  following 
manner 


on  the  first  interior  row  and 
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figure  26.  Vertical  structure  of  the  NORAPS  sigma  surfaces  with  distribution  of 
variables  (Hodur  1987). 


Kh  =  —  1.25V.r  (.1-2) 

on  all  other  interior  points,  where  Vx  is  the  grid  spacing  in  meters.  A  Robert  filter  with 
a  smoothing  coefficient  of  0.15  is  then  applied  to  control  the  high  frequency  time  oscil¬ 
lations. 

The  NORAPS  model  uses  one-way  influence  boundary  conditions  that  specify  the 
time-dependent  lateral  boundary  values  on  the  finer  mesh  NORAPS  model  from  the 
NOGAPS  predictions.  Although  the  N'OGAPS  solution  forces  the  fine-mesh  model,  the 
fine  grid  has  no  effect  on  the  coarse-grid  solution.  Because  of  timeliness  restrictions  due 
to  operational  restraints,  the  boundary  conditions  are  acquired  from  the  previous  12-h 
forecast.  A  scheme  developed  by  Pcrkey  and  Kreitzberg  (1976)  is  used  to  spatially  in¬ 
terpolate  the  solutions  near  the  boundaries  of  the  NORAPS  model.  The  NOGAPS  time 
tendencies  are  combined  with  the  NORAPS  time  tendencies  over  a  distance  of  several 
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Figure  27.  Horizontal  structure  of  the  NOR APS  with  distribution  of  variables 
(Hodur  1987). 

grid  points  to  reduce  false  reflections  at  the  regional  model  boundary  due  to  the  change 
in  grid  spacing. 

The  numerical  model  would  not  be  complete  without  some  type  of  numerical  simu¬ 
lation  or  the  planetary  boundary  layer  (PBL).  The  PBL  is  defined  as  the  lowest  layer  in 
the  model  atmosphere  and  is  well  mixed  in  temperature,  momentum  and  moisture. 
Interactions  that  occur  between  the  I’BL  and  the  upper  atmosphere  provide  sources  and 
sinks  for  momentum,  heat  and  moisture.  According  to  Hodur  (1988)  NORAPS  makes 
use  of  a  multilayer  PBL  scheme  that  is  based  on  Louis  (1979)  and  is  a  change  from  the 
PBL  scheme  described  by  Hodur  (1987). 

According  to  Hodur  ( 1987),  the  NORAPS  uses  a  modified  version  of  the  Kuo  (1965) 
cumulus  parameterization  method.  This  scheme  links  the  convection  to  PBL  variables 
by  requiring  moisture  convergence  in  the  PBL.  The  moisture  convergence  is  defined  as 


(A-  3) 


M,  =  y  V qmnvm  x  (1  -  <rPBL)  +  pswq\, 

where  the  first  term  on  the  right  side  is  the  vertically-integrated  moisture  convergence 
and  the  second  term  is  the  surface  moisture  flux.  Convection  occurs  when 
M,  >  7.055  x  and  when  the  equivalent  potential  temperature  decreases  with 

height  from  the  PBL  to  the  first  layer  above  the  PBL.  The  final  constraint  to  this 
process  is  that  the  convection  cannot  occur  if  the  lifting  condensation  level  is  above  the 
PBL,  which  insures  that  the  mixing  is  strong  enough  to  force  convection.  This  scheme 
separates  the  moisture  transport  into  two  parts,  as  described  by  Kuo  (1974).  A  value 
of  b  is  defined  such  that  (bM,)  is  used  to  moisten  the  environment  and  (1  —  b)M,  con¬ 
denses  and  falls  instantaneously  as  rain.  The  factor  b  is  defined  as  the  vertical  average 
of  1  —  relative  humidity  (RH),  where  RH  is  defined  between  0  and  1.  Introduction  of 
the  factor  b  tends  to  moisten  the  column  more  for  a  low  RH  and  to  condense  more 
moisture  for  a  high  RH.  The  modified  Kuo  scheme  provides  only  the  physics  of 
convective  precipitation,  but  does  not  describe  large  scale  (non-convective)  precipitation. 

N on-convective  precipitation  is  assumed  to  occur  when  supersaturation  conditions 
occur  at  a  gridpoint.  Moisture  is  condensed  until  saturation  is  reached,  and  the  con¬ 
densed  moisture  falls  into  the  next  layer  below  and  will  evaporate  until  saturation  oc¬ 
curs.  If  the  layer  below  is  already  super-saturated,  the  moisture  will  continue  to  fall  to 
the  next  layer.  Thus,  precipitation  will  occur  only  if  the  air  is  saturated  from  the  cloud 
to  the  ground.  The  precipitation  routines  are  called  every  eight  time  steps,  which  en¬ 
hances  computer  efficiency.  The  heating  and  moistening  rates  are  then  spread  evenly 
over  the  subsequent  eight  time  steps. 

The  NORAPS  forecast  model  also  takes  into  account  solar  radiation,  which  is  es¬ 
sential  in  the  prediction  of  surface  temperatures,  and  calculates  cooling  rates  of  cloud 
tops  that  may  deepen  cloud  layers.  The  radiation  parameterization  in  the  NORAPS 
model  follows  Katayama  (1974)  for  short-wave  radiation  and  Sasamori  (1968)  for  long¬ 
wave  radiation. 

The  final  component  of  NORAPS  is  the  output  of  the  forecast  fields  for  use  by  the 
forecaster.  Data  can  be  interpolated  to  1 5  pressure  levels.  Recent  changes  in  the  model 
have  improved  the  resolution  below  850  mb  and  above  100  mb.  Output  parameters  in¬ 
clude  winds,  temperature,  surface  pressure,  relative  humidity,  absolute  vorticity  and  la¬ 
tent  heat  release.  The  forecaster  has  the  option  of  using  Mercator,  Lambert  conformal 
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or  stereographic  projections,  which  helps  minimize  distortions  in  the  tropics, 
midlatitudes  or  polar  regions,  respectively. 

A  bogus  vortex  has  been  incorporated  into  the  NORAPS  due  to  the  relative  lack 
of  data  in  the  vicinity  of  tropical  cyclones  (Hodur  1988).  The  bogus  vortex  is  derived 
from  a  model-generated  storm  produced  by  integrating  NORAPS  to  a  steady-state  from 
an  initially  symmetric  circulation,  with  a  constant  sea-surface  temperature  and  a  no¬ 
mean  flow  environment.  The  gridpoints  from  this  forecast  are  used  as  observations  in 
a  successive-correction  analysis  to  incorporate  the  storm  structure  into  the  initial  fields. 
These  model  observations  are  dominant  within  800  km  of  the  center  and  are  linearly 
reduced  to  zero  at  1600  km.  A  heirarchy  of  storms  have  been  generated  based  on  the 
initial  latitude.  The  spin-up  storm  that  most  closely  matches  the  latitude  of  the  storm 
is  selected  for  the  bogus.  Prior  to  the  1988  season,  a  change  was  made  to  use  the  beta- 
plane  rather  than  the  f-plane  in  the  model-generated  storm  (Hodur  1989). 
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Table  1.  NORAPS  MODEL  SIGMA  LEVELS. 


Level  number 

Original  sigma  values 

Intermediate  sigma  values 

1 

.02 

.01 

2 

.05 

.035 

3 

.085 

.675 

4 

.125 

.105 

5 

.175 

.15 

6 

.225 

.20 

7 

.275 

.25 

8 

.325 

.30 

9 

.375 

.35 

10 

.450 

.4125 

11 

.550 

.50 

12 

.650 

.60 

13 

.750 

.70 

14 

.825 

.7875 

15 

.875 

.85 

16 

.905 

.89 

17 

.935 

.92 

18 

.965 

.950 

19 

.985 

.975 

20 

.995 

.990 

21 

1.00 

.9975 

\ 
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APPENDIX  B.  DATA  ACQUISITION  AND  PROCESSING 

A.  DATA  ACQUISITION 

The  NORAPS  analysis  and  forecast  fields  were  obtained  from  Dr.  Rich  Hodur  of 
the  Naval  Environmental  Prediction  and  Research  Facility  (NEPRF).  Data  from  the 
9-track  tape  were  transferred  to  a  mass  storage  device  on  the  Naval  Postgraduate  School 
(NPS)  IBM  3033  computer.  The  data  fields  are  on  a  79  x  69  staggered  C-grid  with  21 
sigma  levels.  Sigma  is  the  vertical  coordinate  used  by  the  model  and  is  defined  as 

p 

a  =  — ,  where  p  is  the  pressure  at  any  height  and  ps  is  the  terrain  pressure.  Sigma  is 
defined  as  1  at  the  surface  where  p  =  ps  and  zero  at  the  top  of  the  model  where  p  —  0. 
The  horizontal  grid  spacing  used  in  the  model  is  80  km.  The  analysis  time  for  the  model 
forecast  is  0000  UTC  4  October  1988.  A  72-h  forecast  is  generated  with  output  fields 
available  every  12  h  until  0000  UTC  7  October  1988. 

The  NORAPS  model  output  fields  are  on  a  Lambert  conformal  projection.  This 
scheme  is  a  bi-conic,  secant  type  of  projection  that  preserves  angles  when  projecting  the 
earth's  surface  onto  a  plane  surface.  The  Lambert  conformal  projection  scheme  is  best 
suited  for  the  midlatitudes  due  to  the  minimal  distortion  between  the  true  parallels  of 
30°  and  60°N.  When  used  in  the  tropics,  this  scheme  produces  distortion  of  land  masses 
and  representations  of  meteorological  variables  near  the  equator.  These  deficiencies  are 
alleviated  by  modifying  the  graphics  display  program  used  in  this  study. 

The  variables  from  the  model  forecast  are  the  u  and  v  wind  components,  temper¬ 
ature.  specific  humidity,  terrain  pressure,  mean  sea-level  pressure  and  ground  temper¬ 
ature.  Along  with  these  variables,  the  forecast  fields  also  include  the  surface  sensible 
heat  flux,  surface  latent  heat  flux,  convective  and  stable  precipiatation. 

Manipulations  of  the  variables  were  performed  on  the  NPS  IBM  mainframe  com¬ 
puter.  The  results  were  then  transferred  to  main  storage  device  on  the  NPS 
Meteorology/Oceanography  Interactive  Digital  Environmental  Analysis  (IDEA)  Labo¬ 
ratory.  It  was  decided  to  make  use  of  the  IDEA  lab  due  to  the  availability  of  display 
programs  designed  by  Prof.  Wendell  Nuss,  Dept,  of  Meteorology. 

B.  DATA  PROCESSING 

The  concept  of  isentropic  potential  vorticity  is  based  on  the  conservation  of  poten¬ 
tial  vorticity  for  air  flow  on  a  theta  surface  (Hoskins  ei  al.  1985).  Potential  vorticity  may 
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change  due  to  latitudinal  displacements  (Coriolis  parameter)  or  by  diabatically  changing 
the  separation  of  the  isentropic  layers. 

Potential  vorticity  is  computed  on  a  three-dimensional  array  of  grid  points  using  the 
expression  from  Hoskins  et  al.  (1985) 


;>>k=?x  (k,+/i 


dv  86  du  86  \ 
dp  dx  dp  dy  j 


(5-1) 


where  IPV  is  the  potential  vorticity  on  an  isentropic  surface,  is  the  relative  vorticity 
on  a  constant  pressure  surface,  f  is  the  Coriolis  parameter,  6  is  the  potential  temperature, 
p  is  the  pressure  and  u  and  v  are  the  zonal  and  meridional  wind  components.  The  last 
two  terms  on  the  right  side  of  the  equation  are  necessary  for  the  coordinate  transfor¬ 
mation  from  pressure  to  isentropic  coordinates.  Because  the  NORAPS  output  uses 
sigma  as  the  vertical  coordinate,  another  coordinate  transformation  must  be  incorpo¬ 
rated  into  (5  —  1).  The  transformation  from  pressure  to  sigma  coordinates  was  com¬ 
puted  using  the  following  identities  as  described  by  Phillips  (1957), 


dv  1  dv 

dp  Ps  da  ' 

(5-2) 

du  1  du  j 

_________  ancj 

cp  ts  ca 

(5-3) 

do  i  ee 

dp  Ps  da 

(5-4) 

Using  these  identities,  IPV  was  calculated  on  an  isentropic  surface  from  the  sigma  sur¬ 
face  fields. 

The  NORAPS  model  output  is  on  a  staggered  C-grid.  The  major  disadvantage  of 
this  scheme  for  calculating  IPV  is  that  the  u  and  v  wind  components  are  not  located  at 
the  same  gridpoint.  An  averaging  technique  is  used  to  define  values  of  u  (v)  to  the 
north-south  (east-west)  of  a  central  gridpoint  where  6  is  defined 


UMEA\(I,  J)  = 


u(I  —  1,  J)  4-  u(l,  J)  +  u(I  —  1,  J  4-  1)  +  u(I,  J+  1) 

4 


(5-5) 


which  gives  a  mean  value  of  u  at  the  original  v  gridpoints.  A  similiar  procedure  is  used 

«■> 

to  define  the  mean  v  component  for  the  east-west  derivative.  The  values  of  -rr-  and 

dv  8y 

-r-  are  then  used  in 
c.x 
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(B-  6) 


r  & v  du 
p  8x  dy 


to  compute  the  relative  vorticity.  The  values  of  relative  vorticity  at  each  sigma  level  are 

then  vertically  averaged  to  give  a  value  of  (,  at  a  point  between  two  successive  sigma 

levels.  This  is  required  because  -4r~ ,  -rr-  and  are  calculated  between  sigma  surfaces. 

op  op  op  e 

Therefore,  IPV  is  also  calculated  at  an  intermediate  level  between  two  sigma  levels. 

Six  theta  surfaces  are  examined  to  describe  the  key  features  for  this  study.  A  major 
consideration  is  to  assure  that  the  vertical  resolution  is  fine  enough  so  that  even  small 
values  of  IPV  can  be  accurately  determined.  By  using  21  sigma  levels  and  then  trans¬ 
forming  to  theta  surfaces,  a  more  detailed  representation  of  IPV  is  achieved.  This  is  a 
major  advantage  over  earlier  IPV  studies  in  which  mandatory  pressure  levels  were  used 
(Elsberry  and  Kirchoffer  1988).  Detailed  structures  may  be  lost  due  to  the  decrease  in 
vertical  resolution. 


A  linear  interpolation  scheme  is  used  to  transform  the  coordinate  system  from  sigma 
to  theta  surfaces.  A  theta  search  is  defined  with  values  ranging  from  300  to  360  K.  The 
theta  search  value  is  bracketed  between  two  sigma  levels  at  each  grid  point.  The  u  and 
v  wind  components,  pressure  and  IPV  are  then  interpolated  linearly  at  the  defined  sigma 
surface.  The  value  of  these  variables  corresponds  to  a  specific  theta  surface  at  each  grid 
point.  Using  the  IDEA  Lab  and  the  programs  designed  by  Prof.  Xuss,  these  variables 
are  plotted  on  horizontal  theta  surfaces.  Vertical  cross-sections  also  are  used  to  illus¬ 
trate  the  vertical  structure  of  these  variables  in  the  model  analysis  and  forecast. 


V 
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